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Abstract 

In this work, we have carried out the influence of temperature dependent viscosity on thin film flow of a 
magnetohydrodynamic (MHD) third grade fluid past a vertical belt. The governing coupled non-linear differential 
equations with appropriate boundary conditions are solved analytically by using Adomian Decomposition Method (ADM). 
In order to make comparison, the governing problem has also been solved by using Optimal Homotopy Asymptotic 
Method (OHAM). The physical characteristics of the problem have been well discussed in graphs for several parameter of 
interest. 
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Introduction 

The subject of non-Newtonian fluids is popular and is an 
area of active research specially in mathematics, industry and 
engineering problems. Examples of non-Newtonian fluids 
include plastic manufacturing, performance of lubricants, food 
processing, movement of biological fluids, wire and fiber 
coating, paper production, transpiration cooling, gaseous diffu- 
sion, drilling mud, heat pipes etc. These fluids are described by 
a non-linear relationship between stress and the rate of 
deformation tensors and therefore several models have been 
proposed. There are several subclasses of non-Newtonian fluids. 
Third grade fluid is one of the important fluid in this category 
and its equation is based on strong theoretical foundations, 
where relation between stress and strain is not linear. Therefore, 
in this problem, we have considered third grade fluid. 
Considerable efforts have been made to study non-Newtonian 
fluids for various geometrical configurations via analytical 
techniques. Some developments in this direction are discussed 
in [1—19]. On the other hand, the physical importance of thin 
film has been highlighted by scientists and engineers. Amongst 
them, Khalid and Vafai [20] studied hydrodynomic squeezed 
flow and heat transfer over a sensor surface. Miladinova et al. 
[21] investigated thin film flow of a power law liquid falling 
from an inclined plate where it was observed that saturation of 
non-linear interaction occur in a permanent finite amplitude 
wave. 

Similarly, Taza Gul et al. [22] investigated effects of slip 
condition on thin film flow of third grade fluids for lifting and 
drainage problem under the condition of constant viscosity. The 
effects of various parameters on the lift and drainage velocity 
profiles are also studied. 



It is crystal clear that the physical problems are frequently 
modeled, using non-linear differential equations. Recently, several 
analytical and numerical techniques were used for solution of such 
non-linear problems. In order to find analytical approximate 
solutions of non-linear differential equations, researchers usually 
use approximate techniques such as Homotopy Perturbation 
Method (HPM) [23], Homotopy Analysis Method (HAM) [24] 
and Optimal Homotopy Asymptotic Method (OHAM) [25]. 
OHAM is a powerful mathematical technique and has already 
been applied to several non-linear problems. Marinca and 
Herisanu [26] used OHAM for solving non-linear equations 
arising in heat transfer problems. In another paper, Marinca [27] 
applied OHAM to study steady flow of a fourth grade fluid past a 
porous plate. Joneidi et al. [28] analyzed micropolar flow in a 
porous channel with high mass transfer. Siddiqui et al. [29] 
examined a thin film flow of non-Newtonian fluid over a moving 
belt. In another study, Siddiqui et al. [30] discussed the thin film 
flow of a fourth grade fluid down a vertical cylinder. Costa and 
Macedonio [31] noticed that increase in velocity may produce 
additional growth of local temperature. Nadeem and Awais [32] 
investigated thin film unsteady flow with variable viscosity. They 
analyzed the effect of variable thermo capillarity on the flow and 
heat transfer. Ellahi and Riaz [33] discussed analytical solution for 
MHD flow in a third grade fluid with variable viscosity. Whereas 
Aksoy et al. [34] found an approximate analytical solution for flow 
of a third grade fluid through a parallel plate channel filled with a 
porous medium. 

The main objective of this research is to study thin film flow of 
MHD third grade fluid over a vertical belt under the influence of 
temperature with variable viscosity. More exactly, we are 
interested in showing the effects of MHD and variable viscosity 
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Figure 1. Geometry of the problem (a) Lift problem and (b) Drainage problem. 

doi:10.1371/journal.pone.0097552.g001 



with heat transfer in a thin film fluid flow such as silicate melts and 
polymers. In these fluids, viscous friction generates a local increase 
in temperature near the belt with decrease in resultant viscosity 
and frequently increases the flow velocity. The governing problem 
is solved using an analytical technique known as Adomian 
Decomposition Method (ADM). This technique was introduced 
by Adomian [35,36] for finding the approximate solutions for 
linear and non-linear differential equations. Wazwaz [37,38] used 
ADM for reliable treatment of Bratu-type and Rmden-Fowler 
equations. For comparisons and accuracy of results, the governing 
problem has also been solved by using OHAM. 

Basic Equations 

The continuity, momentum and energy equations for incom- 
pressible, isothermal and electrically conducting third grade fluid 
are; 
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Here, p is the constant density, g denotes gravitational 
acceleration, u is the velocity vector of the fluid, 0 defines 
temperature, K is the thermal conductivity, C p is specific heat, 

L = Vu, — = — + (w.V) denotes material time derivative, / is the 




Figure 2. Comparison of ADM and OHAM methods for lift 
velocity profile. St =0.1,M 2 =0.2,^=0.6,^ = 0.1,5, =0.4, 
Ci =0.212619, Ci = -0.243971, C 3 =0.016110, Ca = -0.0168232. 
doi:1 0.1 371 /journal.pone.0097552.g002 



Figure 3. Comparison of ADM and OHAM methods for lift 
temperature distribution. St =l,Af 2 =0.2,/?= \.2,A = 0.01, B r = 
4, Ci =-0.9489619,0 = -0.000257, Ci =-13.262703, C 4 = 
-163.37897. 

doi:10.1371/journal.pone.0097552.g003 
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Figure 4. Comparison of ADM and OHAM methods for 
drainage velocity profile, s, =0.09, M 2 = 0.1,0 = 0.6,^1 = 0.01, 
B r =0.4, Ci =-0.9489619, Ci = -0.000257, =- 13.262703, C 4 = 
-163.37897. 

doi:1 0.1 371 /journal.pone.0097552.g004 

current density and T is the Cauchy stress tensor. Moreover, a 
uniform magnetic field _B=(0, Bofi), is applied in a direction, 
perpendicular to the belt. The Lorentz force per unit volume is 
given by 



jxB=(0,aBlu(x),0), 
The Cauchy stress tensor T, is given by 



(4) 



Figure 6. Influence of the Brinkman number on the lift velocity 
profile when S, = 1,/1 = 0.1, Af 2 =0.5,0 = 2.0. 
doi:10.1371/journal.pone.0097552.g006 



Here a.,. ,„,,/?, , „ „, are the material constants and An A, A 2 , 
./(' = 1,2) "y(7=l,2,3) u > '• " 

A 3 are the kinematical tensors given by 



Ao =1 

Ai = (V«) + (V« 
DAi 
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Dt 
DAi 
Dt 



Ai(Vu)+(Vuy Ai , 



Ai(Vu)+ (Vu)' Az, 



(7) 



T=- P I + t (5) 

where —pi denotes spherical stress, p is the hydrostatic pressure 
and shear stress tensor x, is defined as 

t = HAi + ai Ai + o( 2 A\ + P\ Ai 

(6) 

+ fa{A\ Ai + AzA\) + fa{trAty A\ , 
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Figure 5. Comparison of ADM and OHAM methods for 
drainage temperature distribution. Si =0.1, M 2 =0.3,0 = 0.6, 
A = 0.6,B, =0.3, Ci =-1.024309,C2 = -2.211328, Ci =-0.002455, 
C 4 =0.000389. 

doi:1 0.1 371 /journal.pone.0097552.g005 



Formulation of Lift Problem 

Consider, a wide flat belt moves vertically upward at a constant 
speed U through a large bath of third grade liquid. The belt carries 
a layer of liquid of constant thickness, S with itself. For analysis, 
coordinate system is chosen in which the y-axis is taken parallel to 
the surface of the belt and x-axis is perpendicular to the belt. 
Uniform magnetic field is applied transversely to the belt. It is 
assumed that the flow is steady and laminar after a small distance 
above the liquid surface layer and the external pressure is 
atmospheric everywhere. 
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Figure 7. Influence of the Brinkman number on the lift 
temperature distribution when Si = \,A = QA, M 2 =0.5,0 = 2.0. 
doi:10.1371/journal.pone.0097552.g007 
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Figure 8. Effect of Stock number on the lift velocity profile, 
when M 1 = 0.2,0 = 1.2,^1 = 0.01, ft =4.0. 
doi:10.1371/journal.pone.0097552.g008 



Figure 10. Effect of viscosity parameter on the lift velocity 
profile, when M 2 =0.2,5, = 0.4,,8 = 0.6, ft =4.0. 
doi:10.1371/journal.pone.0097552.g010 



Velocity and temperature fields are 

k=(0,w(x),0) and © = ©(x) 



(8) 



(12) 



Using the velocity field given in Eq. (8) the continuity Eq. (1) 
satisfies identically and Eq. (5) gives the following components of 
stress tensor: 



Txz = Tyz =0, 



(13) 



Txx = -p + (2ai + cc 2 ) 



du 

~dx) ' 



Incorporating Eqs. (9-13) into the momentum and energy 
(9) equations (2, 3), we get 
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Figure 9. Effect of Stock number on the lift temperature 
distribution when M 2 =0.2,/?= l.2,A = 0.01, ft =4.0. 
doi:1 0.1 371 /journal.pone.0097552.g009 



Figure 1 1 . Effect of viscosity parameter on the lift temperature 
distribution when M 2 =0.2, Si =1,/?= 1.2, ft =50. 
doi:10.1371/joumal.pone.0097552.g011 
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Figure 1 2. The effect of magnetic force on lift velocity profile 
when ft. =4.0,4 = 0.1,0 =2.0, St =0.4 
doi:10.1371/journal.pone.0097552.g012 



Figure 14. The effect of Stock number on velocity for drainage 
problem when M 1 =0.2,(8 = 1.2,4 = 0.1, ft =20.0. 
doi:10.1371/journal.pone.0097552.g014 



The corresponding boundary conditions are: 
u=U,®= 0 O at.Y = 0, 



(16) 



For Reynold's model, the dimensionless viscosity 



(19) 



du 
dx 



= 0,0= 0i atx = <5, 



Introducing the following non-dimensional variables 
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Using Taylor series expansion, one may represent viscosity and 
its derivative as follows: 
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where B, is the Brinkman number , M 2 is the magnetic parameter, 
p is the non-Newtonian parameter and S t is the Stock's number. 



Using the above dimensionless variables into Eqs. (14-17) and 
dropping out the bar notations, we obtain. 
Eq. (14) has been rectified for Eq. (21, 22) 
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Figure 13. The effect of magnetic force on lift temperature 
distribution when ft =4.0,A =0.1,(8=2.0, St =0.4. 
doi:10.1371/journal.pone.0097552.g013 



Figure 15. The effect of Stock number on temperature for 
drainage problem when M 2 =0.2,(8= 1.2,4 = 0.1, ft =20.0. 
doi:10.1371/journal.pone.0097552.g015 
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I® 

dx 2 



+ Br 



du\ 2 ^ / du\ 2 ^ / du s 4 
(ix/ \dx) \dx 



w=l,0 = O at x = 0, 



du 

dx 



= 0,0 = 1 at x = 0, 



= 0, (22) 

(23) 
(24) 



Solution of Lifting Problem 

The OHAM solution. In order to solve the system of 
equations (21-24), we define the linear, non-linear functions and 
source terms respectively as follows: 



LlW 



dx 2 
+ 6P 



,Ni(x) 



( du d® 



, d u 



\dx dx dx 2 



du\ 2 d 2 u 



(25) 



dx J dx 2 



M 2 u, Gi (x) -- 



S f , 




L2(x)=^,Ni(x) 



= Br 



\dx I \dx 



,&(*)=a 



(26) 



Now, OHAM is applied to non-linear coupled ordinary 
differential Eqs. (21, 22) and Eqs. (25, 26) as follows: 

[1 -P\ [Li (x) + Gi (x)] - Hi {p) [Li {x) + Ni (x) + Gi (x)] = 0, 



[ 1 -p] [L 2 (x) + Gi (*)] - H 2 (p) [L 2 (x)+N 2 (x) + Gi (x)] = 0, 

We consider u(x),®(x), Hi{p), H 2 {p) as the following. 

Hi 0) =p Ci + P 2 C 2 , H 2 (p) =p C 3 + P 2 Ca ,u{x) 
= u 0 (x) +pu x (x) + P 2 u 2 (x) , 

0(x) = 0 o (x) +p®i{x) + P 2 ® 2 {x), 



ill) 



(28) 



Substituting u(x),®(x), Hi (p), H 2 (p), Li (x), L 2 {x),Nl (x), 
N 2 {x), Gi(x), G 2 {x) from Eq. (28) into Eq. (27) and after some 
simplifications based on power of p-terms, we get the following. 

Zero components: 



P": -S, +^=0,^=0, 
dx 1 dx 1 



(29) 



Figure 17. Drain temperature distribution for various values of 
viscosity parameter when M 2 = 0.2,m = 1.0/3 = 0. 6, B r =50.0. 
doi:1 0.1 371 /journal.pone.0097552.g01 7 
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First components: 
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Solving Eqs. (29, 30) along with boundary conditions (23, 24), 
we get the term solutions as fallows. 
Zero term solution: 



u(x) =Ax + B+ St^r + L l [M 2 u]-6PL~ 
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The series solutions of Eqs. (35, 36): 



(36) 



u 0 (x)=-(2-2xSt + x 2 St),®o(x) = x, (31) 



First term solution: 
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The second term solution for velocity and temperature are too 
bulky, therefore, only graphical representations up to second order 
are given. 

The series solutions of velocity profile and temperature 
distribution are 



J2u n = Ax + B+ S,^ + M 2 L~ l \Y^ u n J -6jS.fr 1 



+ AL~ 
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(38) 



The Adomian polynomials An . B n , C« , Z>« and _F„ , for Eqs. 
(37, 38) are defined as 



n = 0 



rf« \ 2 fd 2 U 



dx) \dx 2 )' 



(39) 



«(x, Ci ) = «o to + u\ (x) + u 2 (x), 
and ®(x, Ci) = ©o to + 0i to + ©2 to 



(33) 



du\ ( d® 
dx I \dx 



(40) 



The arbitrary constants Ci,i= 1,2,3,4 are found out by using 
the residual 



i? = L[w(x, C;)] + N[u(x, c t )] + G[u(x, c( 



(34) 



(41) 



For velocity profile and temperature distribution the arbitrary 
constants are mentioned in graphs. 

The constants Cl, C2 C3 C4 can also be obtained from 
Collocation and Ritz methods. 

The ADM solution. The inverse operator = § du, of 
the ADM on the second order coupled Eqs. (21,22) is used: 



(42) 
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The components of Adomian polynomials are derived from 
Eqs. (39-44) as: 



0i(x) = — B,L~ X [Do] + Br A IT 1 [Eo] ~ 2 B, P L~ 1 [Fo] , (52) 
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(46) 



The series solutions of Eqs. (37, 38) are derived as: 

x 2 

U 0 + U\ + U 2 +...=Ax + B+ S,y + 

M 2 L' 1 {u 0 + Ul + u 2 + ...) -6P L'^Ao + A\ + Ai +■■■]+ ^ 
AZT'ltfo + Bi + Bi + ...]+AZT 1 [Co + Ci + C2 +••■], 



©2 (x) = - B r IT 1 [Di] + B r A ZT 1 [Ei] - 2 B r jS IT 1 [Fi] , (54) 
subject to the boundary conditions 



«b(0) = 1, ^ ( 1 ) = 0, «„(0) = 0, ( 1 ) = 0,« > 1 , (55) 



®o(O)=O,0 o (l) = l,© n (O)=O,®„(l)=O,«>l, (56) 



Using boundary conditions from Eqs. (55, 56) into Eqs. (49-54), 
we obtain 



w 0 (x) = -(2-2SVx+svx 2 ), 



0 o (x)=x, 



(57) 



(58) 



00 + 0] +0 2 +... = Ex+F- BrL~ l [Do +Di + D2+-] 

+ RAr'[£o + £i+£2+-] (48) 
-2B r PL- l [Fo+ Fx +F 2 +•••], 



The velocity and temperature components are obtained by 
comparing both sides of Eqs. (47, 48): 

Components of the lift problem up to second order are: 



u 0 (x)=Ax + B+ Sty, 



(49) 



«,(.vj=( -M 1 + K *^+*PsVjx 
+ (^ + ^-6psAx 2 



+ 



+ 



M 2 St , A S, , ,„ r . : . 1 ; 



6 • f -4/;.v;) v 

^-4/^V, 



(59) 
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(11 2 
^B r S 2 , - ^Afi, S 2 + ^PBrSt )x 



+ QaS? +^AB r S 2 +^[tB r st^Jx i (60) 
+ ( - ^ B r S 2 - ^ A B r S 2 , + 2p B r S^) x 4 

+ (^ABrS 2 +^pB r S?)x 5 +(~^PBrSt)x 6 , 



Due to lengthy calculations, the analytical results have been 
given up to first order but they have been shown graphically up to 
second order. 

Formulation of Drainage Problem 

Under the same assumptions as in the previous problem, we 
consider a film of non-Newtonian liquid draining down the 
vertical belt. The belt is stationary and the fluid drains down the 
belt due to gravity. The gravity in this case is opposite to the 
previous case. The coordinate system is selected same as in the 
previous case. Assuming that the flow is steady and laminar, 
external pressure is neglected whereas the fluid shear forces keep 
gravity balanced and the thickness of the film remains constant. 

Boundary conditions for the drainage problem are 



du 

u = 0, at x = 0, -r- =0, at x = <5, 

dx 



(61) 



Using non-dimensional variables, the boundary conditions for 
drainage problem become 



w„(0)=0, and ^0)=0, «>0, 



(62) 



OHAM is applied to non-linear coupled ordinary differential 
Eqs. (63, 64) as 

[1 -p] [Li (x) + Gi M] - Hi (p) [Li (x) + Ni (x) + Gi (x)] = 0, 



[\-p][L 2 {x)+G 1 {x)]- Hi(p)[Li(x)+ N 2 (x) + Gi(x)]=0, (65) 



We consider u(x),®(x), Hi (p), H 2 {p) as the following 

Hi(p)=pCi +P 2 C 2 , 
H 2 (p)=pCi +P 2 C 4 , 
u(x)= u 0 (x) +p U\(x) + P 2 u 2 (x), 
0(x)=0 o (x) +p@i(x) + P 2 ® 2 (x), 



(66) 



Substituting w(x),0(x), Hi (p), H 2 (p), Li (x), Ni (x), Gi(p), 
L 2 {x),N 2 {x) and G 2 (x) from Eq. (66) into Eq. (65) we have the 
following components of velocity and temperature. 

Zero components: 



dx 2 



dx 2 



(67) 



First components: 



/»>: -S.-S.C^mAaCi^^-^-Ci^ 
ax dx dx dx L dx A 



+ ACi ® 0 (x] 



d 2 »o 
dx 2 



-npCi 



! 'no\ 2 d 2 uo d 2 ui 



= 0, 



_ BrcJ^\ 2 +ABrCnol im its 3 e 0 (x) '-Aft W^Y 



~ 4[i ^\-dx-J 



d 2 ®o n d 2 ®o , d 2 ®i 



dx 2 



+ - 



dx 2 dx 



= 0, 



(68) 



For temperature distribution, the boundary conditions are same 
as given in Eq. (56). 

Solution of Drainage Problem 

The OHAM solution. From Eqs. (21, 22), the linear, non- 
linear functions and source term (in drainage case), it is opposite 
due to gravity) are respectively defined as 



d 2 u 



Li = " " Ni = -\(— — + 0—^+66 (^] 2 — 

dx 2 ' \dx dx dx 2 ) \dxj dx 2 (63) 



- M u, Gi = St, 



Li = 



d 2 u 
dx 2 ' 



v.-*l(*)-Ae(±) +V (i)W,-.> 



(64) 



Solving Eqs. (67, 68) with boundary conditions (61, 62), we get 
u 0 (x)=-(2xSt - x 2 S,),®o(x) = x, (69) 



ui(x) 



Ci 
24 



(8M 2 5, + 96/?,S' 3 )x+(144/? > S l3 -12AS,)x 2 
+ {^AS,~4M 2 S 1 +96l3sf)x 3 
+ {M 2 S, -24pS?)x 4 



®i(x) (70) 
■ (3 Sj A- 15 Sj -40/? Sf)x+ (30 S 2 + 120 j8 S 4 ) x 2 " 

— (20 sj + 160/? S 4 +10S 2 A) x 3 
+ (5 Sj + 120/3 S 4 + 10 sj A) x 4 

- (48/? S 4 + 3 Sj A) x 5 + (8/J Sj) x 6 



BrCj 

60 
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Like previous problem, results up to first order terms haven 
been obtained. 

ADM solution. Using ADM on Eqs. (21, 22), the Adomian 
polynomials in equations (45, 46) for both problems are same 
whereas the different velocity components are obtained as: 

Components of the Problem 

The boundary conditions of first and second components for 
drainage velocity profile are same as given in Eq. (56). Also, the 
boundary conditions for temperature distribution are same as 
given in Eq. (57) but solution of these components is different, 
depends on the different velocity profile of drainage and lift 
problems. Due to lengthy analytical calculation, solutions up to 
first order terms are included whereas the graphical representa- 
tions up to second order terms are given. Using boundary 
conditions (62) and (56) into Eqs. (49-54), the components, 
solution are obtained as: 

u 0 (x)= X -{2xS, -x 2 S,), (71) 



©o (*)=*, (72) 
u x {x)=(- X -M 2 S, -4psnx+ (^+6/?S, 3 ) x 2 

(73) 

©i (*) = Q Br S 2 -^A Br S 2 + l 0B, sA x 
+ (-~BrS 2 -2j3B r S*)x 2 

+ 0 Br Si+^ABrSi+l pB r <A a- 3 (74) 
+ (-^BrSf-^AB r Sf-2 jW r S 4 ) X 4 

Results and Discussion 

The effect of Stock number St magnetic parameter M 2 , 
Brinkman number Br non-Newtonian parameter /J and viscosity 
parameter A in lifting and drainage problems together with the 
physical interpretation of the problem have been discussed in 
Figs. 1-20. Fig. 1 shows the geometry of lift and drainage 




0.0 0.2 0.4 0.6 0.8 1.0 



x 

Figure 19. The influence of non-Newtonian parameter />' on 
temperature for drainage problem when 

M 2 =0.5,5/ =M = 0.01,.B, =50.0. 
doi:1 0.1 371/joumal.pone.0097552.g01 9 

problems. A comparison of the ADM and OHAM solutions is 
shown in Figs. 2-5 for various values of physical parameters. It is 
found from these figures that ADM and OHAM solutions are in 
good agreement. Figs. 6 and 7 provide variation of velocity and 
temperature distribution for different values of Brinkman number. 
It has been found that velocity decreases whereas temperature 
inside the fluid increases by increasing B, while keeping the other 
parameters fixed. In Fig. 8, we observed that velocity decreases 
with an increase in the Stock number St- Physically, it is true as 
increasing Stock number causes the fluid s thickness and reduces 
its flow. The effect of Stock number St on temperature distribution 
has been illustrated in Fig. 9. It is observed that temperature 0 
increases monotonically for large values of Stock number St- The 
effect of viscosity parameter A on lift velocity u is shown in Fig. 10. 
It is observed that the speed of flow decreases by increasing A. The 
speed of flow is actually caused by shear s thickening and thinning 
effects due to increase and decrease in viscosity parameter. A 
similar situation is observed in Fig. 1 1 where an increase in 
viscosity parameter A decreases temperature distribution. Here, 
the velocity profiles are parabolic in nature and their amplitudes 
depend on the magnitude of the viscosity parameter A. Variations 
of the magnetic parameter M 2 on lift velocity have been studied in 
Fig. 12. Here, it is clear that the boundary layer thickness is 
reciprocal to the transverse magnetic field and velocity decreases 




0.0 0.2 0.4 0.6 0.8 1.0 



x 

Figure 20. Comparison of the present results with published 
work [22] when S, =0.5, M 2 =0.1,/J = 0.6,/I = 0,a= 1.0. 
doi:1 0.1 371 /journal.pone.0097552.g020 
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Table 1. Numerical results show the comparison of present results with published work [22] when 

St =0.5, M 2 =0.1,j8=0.6,yf=0,a=1.0. 



X 


Present work 


Published work[22] 


Absolute error 


0 


1 


1 


0 


0.1 


0.95336886435638888 


0.95336886432889 


1.11022x10 16 


0.2 


0 Q1 163368 

U.7 1 IUJJUO 


0 91 1 63368 

U. J 1 1 UJJUO 


1 1 1023 *10~ 16 

1 . 1 1 UiJ A 1 U 


0.3 


0.8742813014374999 


0.8742813014375 


0 


0.4 


0.84116332088888 


0.84116332088888 


0 


0.5 


0.8123818359375 


0.8123818359375 


0 


0.6 


0.788190612 


0.788190612 


0 


0.7 


0.76891138782638 


0.76891138782638 


0 


0.8 


0.754865408 


0.754865408 


0 


0.9 


0.74632017043749 


0.7463201 704375 


0 


1.0 


0.7434513888888 


0.7434513888888 


0 



doi:1 0.1 371 /journal.pone.0097552.t001 



as flow progresses towards the surface of the fluid. On the other 
hand, temperature profile as shown in Fig. 1 3 indicates that fluid 
temperature increases with magnetic parameter. Fig. 14 shows 
that velocity increases in drainage flow when Stock number St 
increases. Physically, it is due to friction which seems smaller near 
the belt and higher at the surface of the fluid. Further, it is found 
from Fig. 15 that temperature profile also increases when Si is 
increased. Fig. 16 illustrates the effect of variable viscosity 
parameter A on the drain flow. It is observed that at higher 
values of viscosity parameter A, velocity of the fluid increases 
gradually towards the surface of the fluid. However, it is found 
from Fig. 17, that an increase in viscosity parameter A causes 
gradual decrease in temperature field. The effects of non- 
Newtonian parameter /? on drain velocity have been studied in 
Fig. 18. We observed that an increase in ft raises drain velocity 
profile and decreases temperature profile as shown in Fig. 19. 
Finally for the accuracy purpose the present results are compared 
with published work in [22] in Fig. 20 and in table 1. An excellent 
agreement is found. 



3) Third grade fluid on vertical belt wit surface topography. 

4) Third grade fluid on vertical rotating disc with surface 
topography. 

Conclusion 

In this work, we have investigated the thin film flow non- 
Newtonian third grade fluid due to vertical belt and the fluid was 
subjected to lifting and drainage. Analytical solutions of the lifting 
and drainage problems have been obtained using ADM and 
OHAM. It has been shown graphically that these solutions are 
identical. The results for velocity and temperature have been 
plotted graphically and discussed in detail. It has been observed 
that these solutions are valid not only for small but also for large 
values of the emerging parameters. It has been observed that in 
both cases of lift problem velocity decreases while temperature 
increases with increasing Brinkman number Br- However, in 
drainage problem both velocity and temperature increases. 



Future Work 

We intend to carry out researches in future on third grade fluid 
on vertical belt regarding the following discussions: 

1) Time dependent third grade fluid on vertical belt. 

2) Vogel Model third grade fluid on vertical belt with slip 
boundary conditions. 
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